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ABSTRACT 

We present the first measurement of the spatial clustering of mid-infrared selected obscured and unobscured 
quasars, using a sample in the redshift range 0.7 < z < L8 selected from the 9 deg^ Bootes multiwavelength 
survey. Recently the Spitzer Space Telescope and X-ray observations have revealed large populations of ob- 
scured quasars that have been inferred from models of the X-ray background and supermassive black hole 
evolution. To date, little is known about obscured quasar clustering, which allows us to measure the masses of 
their host dark matter halos and explore their role in the cosmic evolution of black holes and galaxies. In this 
study we use a sample of 806 mid-infrared selected quasars and «250,000 galaxies to calculate the projected 
quasar-galaxy cross-correlation function Wp(R). The observed clustering yields characteristic dark matter halo 

masses of log(Mhaio[/^~'A^0]) = 12.T!^^'1 and 13.3!l;|j4 fo'" unobscured quasars (QSO-ls) and obscured quasars 
(Obs-QSOs), respectively. The results for QSO-ls are in excellent agreement with previous measurements 
for optically-selected quasars, while we conclude that the Obs-QSOs are at least as strongly clustered as the 
QSO-ls. We test for the effects of photometric redshift errors on the optically-faint Obs-QSOs, and find that 
our method yields a robust lower limit on the clustering; photo-z errors may cause us to underestimate the 
clustering amplitude of the Obs-QSOs by at most ~20%. We compare our results to previous studies, and 
speculate on physical implications of stronger clustering for obscured quasars. 
Subject headings: galaxies: active — quasars: general — large-scale structure of universe — surveys 



1. INTRODUCTION 

Supermassive black holes with masses > 1O''M0 are ubiq- 
uitous in_flienucleiofloc^_£aJ[a^^ moderate to high mass 
(e.g., Kormendv & Richstonel[T99^ . It is now well estab- 
lished that most of the total mass in black holes in the nearby 
Universe was accreted in luminous episodes with high Ed- 
dington rates (e.g., Soltan.198 2; Yu & Tremaine 2002) , with 
the growth for massive (Mbh ^ IQ^ Mp,) black holes occur- 
ring predominantly at z > 1 (e.g., iMerloni & Heind 120081 : 
IShankar et al.ll2009l) . These rapidly accreting black holes are 
most readily identified as bright optical quasars with char- 
acteristic broad (> 1000 km s"') emission lines, and lumi- 
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nous continuum emission that can dominate the light from the 
host galaxy, particularly at ultraviolet and optical wavelengths 
(e.g iElvis etal.lll994t IRichards et al.l 120061: ISchneider et all 
l2007h . Optical quasars thus provide powerful tools for trac- 
ing the rapid growth of black holes over cosmic time (e.g., 
iCroom et al.ll2Bo4l IRichards et alMM. iFanet al.ll2006h . 

However, it is increasingly clear that a significant frac- 
tion of the quasar population does not show characteris- 
tic blue continua or broad lines because their nuclear emis- 
sion regions are obscured. Key evidence for the existence 
of obscured (Type 2) quasars comes f rom synthesis models 
of the cosmic X- ray background (e.g., IComastri et al.|[T995t 
iGilh et ani2007al) . as well as direct identification of these ob- 
jects through various observational techniques. These include 
selection of luminous quasars with only narrow op tical lines 
jZakamska et al.ll2003L E00I I200I iRe ves et al.l[2008i) or rel- 
atively weak X-ray emission IPtak et al.l 120061: Vignali et al.l 
I2OO6I 120101) . detection of powerful radio galaxies lacking 
strong nuclear optical contin ua or broad lines (e.g., McCarth^ 
1 1 9931: ISevmour et al.l I2007h. an d detection of X-ray sources 
that are optically fai nt (e.g.,lAlexander et al. 2001; Stern et aT] 
120021: iTreister et al.ll20 04: Mainieri et al. 20051, have hard X- 
ray spectra (e.g.. IVignaU et al.,.2009|,). or have radio bright, 
optica lly weak counterparts (e.g.. iMartmez-Sansigre et al.l 
"200^). 

With the launch of the Spitzer Space Telescope, large 
numbers of obscured quasars can now be efficiently identi- 
fied based on their characteristic (roughly power-law) spec- 
tral energy distributions (SEDs) at mid-infrared (mid-IR) 
wavelengths (Ri3-24 fim). Because mid-IR emission is less 
strongly affected by dust extinction than optical and ultravi- 
olet light, obscured quasars can appear similar to their un- 
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obscured counterparts in the mid-IR, but have optical emis- 
sion characteristic of thei r host galaxies. A number of stud- 
ies usin g mid-IR colors (Lacy et al.'l2004); 'Ster n et al.l 120051: 
iRowan- Robinson et al. 2005; Hickox et al. 200^ hereafter 



HOT), SED fitting ( Alonso-H errero et al . 2006; Donlev et al.1 
120071), or selecting o bjects based o n similarities to mid-IR 
quasar templates (e.g.. lPolletta et al.^ 2006) have been success- 
ful in identifying large numbers of dust-obscured quasars, in- 
dicating that a large fraction, and possibly a majority of rapid 
black hole growth is obscured by dust. 

These large new samples enable detailed statistical stud- 
ies that can explore the role of obscured quasars in galaxy 
and black hole evolution. At present there are a number of 
possible physical scenarios for obscured quasars; in the sim- 
plest "unified models", obscuration is attributed to a broadly 
axisymmetric "torus" of dust that is part of the central en- 
gine, s o obscuration is entirely an orie ntation effect (e.g., 
[Antonucci 1993; Urry & Padovanil ll995). Alternatively, ob- 
scuration may not be due to a central "torus" but to larger 
dust structures such as those predicted during major merg- 
ers of galaxie s (e.g..|Silk & Reeslll998l; ISpringel et ani2005l; 
iHopkins et alJ 12006a). and obscured quasars may represent 
an early evolutionary phase when the growing black hole 
cannot produce a high enough a ccretion luminosity t o ex- 
pel th e surrounding material (e.g.. IHopkins et aDl2008l: Ikind 
120101) . Observations have revealed evidence for obscura- 
tion by a "torus" in some cases and by galactic-scale struc- 
tures in others (e.g., ^Z akamska et al. 2005; Page et al. 2004; 
iMartinez-Sansigre et al.l l2009). and while there are examples 
of obscured quasars that show clear signs of radiative feed- 
back on interstellar gas, it is unclear whether they are driv- 
ing the galaxy-scale outflows invoked in evolutionary models 
jGreene et al.ll201ll) . Thus the physical nature of obscured 
quasars remains poorly understood, and analyses with large 
samples of mid-IR selected quasars will be essential for a 
more complete understanding of rapidly growing, obscured 
black holes. 

One particularly powerful observational tool is spatial 
clustering, which allows us to measure the masses of the 
dark matter halos in which quasars reside. Clustering 
studies of unobscured quasars have shown that the masses 
of quasar host halos are remarkably constant with cosmic 
time, with Mhaio ~ 3 x 1 0^^ Mfr^ over the large red- 
shift range < z < 5 (e.g.,'Porciani et al."2004';'Croom et ajj 
2Q05; Co iletani2007 ; Myers et al. 2007; Shen et al. 20071 
da Angela et all I2008t IPadmanabhan et alj 120091; iRoss et'al] 
20091) . This lack of variation in halo mass implies that the 
bias factor (clustering relative to the underlying dark mat- 
ter) is an increasing function of redshift, since the dark mat- 
ter is more weakly clustered earlier in cosmic time. The 
characteristic Mhaio provides a strong constraint on models 
of quasar fueling by the major mergers of gas-rich galax- 
ies (e.g., [K auffmann & Haehnelt 2000; Spri ngel et al J 120051; 
Hopkins et ah j2006bi). secular instabilities (e.g., iMo et alj 
1998; Bow er et al.ll2006t iGenzel et al.ll2008l) or accretion of 
recycled cold gas from evolved stars ('Ciotti & Os triker|[2007l ; 
ICiotti et al. 2010), and may be related to quasars' role in regu- 
lating star formation and the emergence of the red galaxy pop- 
ulati on in halos of roughly similar mass '-^lO^^-lO'-^ Mp, 
(e.g.. lCoil et"al]|2 008; Brown et al. 2008- |Conrov & Wechslerl 
[20091: iTinker & WetzeLi20 1 0.) . 

Despite the power of clustering measurements in under- 
standing quasar populations, little is known about the clus- 



tering of obscured quasars. Some measurements of lower- 
luminosity AGNs indicate no signific ant difference between 
obscured and unobscured sources (Constantin & Vogele^ 
2006; Li et al. 2006; Gandhi et al. 2006; Mandelbaum et al] 
20091: iGilh et al.. 2009; H ickox et al.l i2009). However, these 



AGNs likely have different physical drivers com pared to pow- 
erful quasars (e.g. . IHopkins & Hemquistl2006l) . For obscured 
quasars at high luminosities (Lboi ~ 10^^ erg s"') and high red- 
shift (z > 1), the clustering has remained largely unexplored. 

In this paper we present the first measurement of the clus- 
tering of mid-IR selected obscured quasars and make direct 
comparisons to their unobscured counterparts. We use a large 
sample of quasars (both obscured and unobscured) in the red- 
shift range 0.7 < z < 1.8 selected on the basis of IRAC col- 
ors by H07, using data from the 9 deg^ Bootes multi wave- 
length survey. We also employ a sample of ?s250,000 galax- 
ies with good estimates of photometric redshift, and measure 
the two-point cross-correlation between q uasars and galaxies . 
We utilize a novel method developed by iMvers etal] (120091 
hereafter M09) to derive the projected real-space projected 
cross-correlation function, making use of the full probability 
distributions for the photometric redshifts. 

Throughout this paper we assume a cosmology with il,^ = 
0.3 and Oa = 0.7. For direct comparison with other works, 
we assume Hq = 70 km s"' Mpc"' (except for comoving dis- 
tances and dark matter halo masses, which are explicitly given 
in terms of h = //o/(100 km s"' Mpc"')). In order to easily 
compare to estimated halo masses in other recent works on 
quasar clustering (e.g.. ICroom et al.ll2005l iMvers et aLll2006l ; 
'da An gela et"ani200llRoss et alJ200 9^. we assume a normal- 
ization for the matter power spectrum of erg = 0.84. Photome- 
try is presented in Vega magnitudes. All quoted uncertainties 
are Icr (68% confidence). 

2. OBSERVATIONS 

The 9 deg^ survey region in Boot es covered by the NOA O 
Deep Wide-Field Survey (NDWFS; ITannuzi & Devil 19991) is 
unique among extragalactic multiwavelength surveys in its 
wide field and uniform coverage using space- and ground- 
based observatories. Extensive optical spectroscopy makes 
this field especially well suited for studying the statistical 
properties of a large number of AGNs (C. Kochanek et al. 
2011, in preparation). Further details of the B ootes data set 
have been presented in previo us papers (e.g., iHickox et al.l 
|2007, 2009; Ashby etalJ|2009D- 

Redshifts for this study come from the AGN and Galaxy 
Evolution Survey (AGES; Kochanek et al. 2011, in prepa- 
ration) whic h used the Hectospec multifiber spectrograph on 
the MMT ( Fabricant et all 120051) . We use AGES Data Re- 
lease 2 (DR 2), which includes all the AGES spectra taken in 
2004-2006. Details of the AGN redshifts ai-e given in H0"7l 
and H ickox et al. (2009). 

Optical photometry from NDWFS was used for the selec- 
tion of AGES targets and to derive optical colors and fluxes 
for AGES sources. NDWFS images were obtained with the 
Mosaic- 1 camera on the 4-m Mayall Telescope at Kitt Peak 
National Observatory, with 50% completeness limits of 26.7, 
25.0, and 24.9 mag, in the Bw, R, and / b ands, respectively. 
Photometry is derived using SExtractor (iBertin & ArnoutsI 
[1996) . 

Mid-infrared observa tions are taken from the Spitzer IRAC 
Shallow Survey (ISS; 'Eisenhardt et al." '2004), and Spitzer 
Deep Wide-Field Survey (SDWFS; Ashbv et al...2009.) . ISS 
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Fig. 1. — (a) IRAC color-color diagram showing the selection of the quasar samples using the criteria of IS tem et al.l ( I2005h . The grayscale shows the density 
of sources detected at > 5o" significance in all four bands in IRAC Shallow Survey data. Blue stars and red circles show the QSO- 1 and Obs-QSO samples, 
respectively. The Stern et al. (2005) color-color selection region is shown by the dashed fine, (b) Illustration of the optical-IR color-selection criteria for dividing 
the IR-selected QSO samp le in to unobscured (QSO-1) and obscured (Obs-QSO) subsamples. Shown is observed R — [4.5] color versus bolometric luminosity, 
calculated as described in j) l3.1l Contours show the distribution for all th e H07 IR-selected quasars, while blue stars and red circles show the QSO- 1 and Obs-QSO 
subsamples at 0.7 < z < 1.8 used in this analysis as described in § 13.11 The right and top panels show histograms of color and Lboi, respectively, for the QSO-Is 
(blue solid line) and Obs-QSOs (red dashed line). The contours and color histograms show that a simple cut in optical-IR color clearly separates the QSO samples 
into two populations, while the L|,oi histograms demonstrate that the two samples are very closely matched in luminosity. 



covers the full AGES field in all four IRAC bands (3.6, 4.5, 
5.8, and 8 /zm), with 5a flux Hmits of 6.4, 8.8, 51 and 50 
/ijy respect ively. The IRAC ph otometry for ISS is described 
in detail in iBrodwin et al.1 (1200 6). The more recent SDWFS 
exposures extend these limits to 3.5, 5.3, 30, and 30 /iJy, re- 
spectively. As discussed below, the quasar sample (as defined 
m Eol was selected using ISS data, while the galaxy sam- 
ple for cross-correlation is selected from the full SDWFS data 
set. In computing bolometric luminosities for the quasars, we 
also make use of 24 /im flux measurements available from the 
Multiband Imaging Photometer for Spitzer (MIPS) GTO ob- 
servations (IRS GTO team, J. Houck (PI), and M. Rieke) of 
the Bootes field. Significant fluxes (> 3a) were obtained for 
97% of the quasars in our sample that lie in the region covered 
by MIPS. 

3. QUASAR AND GALAXY SAMPLES 

Our primary analysis is the two-point cross-correlation be- 
tween mid-IR selected quasars and galaxies. In this section 
we give details of the the quasar (both obscured and unob- 
scured) and galaxy samples. 

3.1. Quasar sample 

The quasar£3 are taken from the sa mple of luminous mid- 
IR selected AGNs presented by IHOTI Quasars are identified 
on the basis of their colors in the mid-IR as o bserved by 
Spitze r IRAC, using the color-color criterion of iStern et all 
(12005") (Figure [Di), and are selected such that their best es- 
timates of redshift are at z > 0.7. To the relatively shallow 
flux limits of the IRAC Shallow Survey, the AGN sample is 
highly complete and suffers little contamination from star- 
for ming galaxies (as discus sed in detail in § 7 of IH07I: see 
also lAssef et al.ll20Toll20TTh . 

We note that while lHOTl refers to the sample as "AGNs", their bolometric 
luminosities are estimated to be in the range lO^^-ltf*^ erg s"', coiTespond- 



IH07I showed that at the ISS flux hmits, the IR-selected 
quasars show a bimodal distribution in optical to mid-IR color. 
The selection boundary at R - [4.5] = 6.1 can be interpreted 
as dividing quasars into unobscured (optically bright and so 
"blue" in R-[4.5]) and obscured (optically faint and so "red" 
in R-[4.5]) subsets (Figure [TJi). For the purposes of this 
study these objects will be referred to as "QSO- Is" and "Obs- 
QSOs", respectively; the reader is reminded that the selection 
is based not on optical spectroscopy but only on optical to 
mid-IR color This selection yields samples of 839 QSO- Is 
and 640 Obs-QSOs at z > 0.7. 

A detailed study of the optical colors, morphol ogies, and 
average X-ray spe ctra o f these objects is given in IH071 To 
briefly summarize, 'H07' found that the QSO- Is have blue op- 
tical colors, point-like optical morphologies, and soft X-ray 
spectra characteristic of unobscured quasars, while the Obs- 
QSOs had redder optical colors, extended optical morpholo- 
gies, hard X-ray spectra and high Lx characteristic of ob- 
scured quasars. The sample does not include all obscured 
quasars, as sources with very large extinction may fall below 
the IR flux limits of the survey or move out of the Ster n et all 
(|2005[) selection re gion (as shown in Figure 1 of hotT see 
also iGoriian eFaiH^OOl lAssef et alll2010h . The typical ab- 
sorbing column for the Obs-QSO sample is estimated to be 
A^H ^ 1022-10^3 cm-2. We expect the Obs-QSOs to suffer 
little contamination from bright star-forming galaxies. IHOTI 
used an X-ray stacking analysis and constraints from deeper 
surveys to estimate the possible contamination, and concluded 
that the contamination is at most «30%, and likely signifi- 
cantiy smaller (<10%). 

i ng roughly to an X-ray luminosity range 5 X 10"*' ^ X lO*' erg s"' 

iMarcoiii et al JI2004t iHopkins et 31.1120071) . Such high luminosity AGNs are 
typically referred to as "quasars" in the literature, so to avoid confusion with 
studies of lower-luminosity active galaxies, here we refer to our sample as 
"quasars". 
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Fig. 2. — Left: Sky positions of IR-selected quasai's at 0.7 < ; < 1.8 in the Bootes field. Quasars ai'e selected using the color-color ciiteria of Stem et al.' J2005|), 
and are divided into unobscured (QSO-1; blue stars) and obscured (Obs-QSO; red circles) as sources with optical-IR color blueward and redward, respectively, 
of the boundary S — [4.5] =6.1. Right: Surface density on the sky for the sample of 256,124 SDWFS galaxies at 0.5 < z < 2. Objects in regions of unreliable 
photometry are excluded from both the quasar and galaxy data sets. 



For our spatial correlation analysis, we limit the IR-selected 
quasar sample to the redshift range 0.7 < z < 1.8, to max- 
imize overlap with the normal galaxies in the field (§ 13. 2t . 
We also include only objects in regions of good optical pho- 
tometry and away from bright stars. These criteria yield 
563 QSO-ls and 361 Obs-QSOs. Finally, we restrict the 
QSO-1 sample to those spectroscopically identified as broad- 
line AGNs, to ensure that they unambiguously represent a 
sample of unobscure d qua sars and to enable clean tests of 
photo-z errors (see § I6.3t . Of the full sample of QSO-ls 
all redshifts, the vast majority (80%) have optical spectra 
from AGES and 96% of these are classified as broad-line 
AGNs at 0.7 < z < 4.3, supporting their selection as un- 
obscured quasars. We limit the QSO-1 sample to the 445 
that have accurate optical spectroscopic redshifts in the range 
0.7 < z < 1.8 and clear broad emission line features. (In 
a sense this is conservative; we verify that including the 
20% of objects with only photo-zs has no significant effect 
on the clustering results.) Based on these selection criteria, 
our QSO-1 sample is essentially equivalent to other Type 1 
quasar samples selected p urely on optical photometric colors 
and/or spectroscopv (e.g.. [ R ichards et al."2001\ ' Croom et al.l 
l2004t [Schneider et alj|2007l; iRichaids et al. 2009157 since the 
vast majority of spectroscopic Type 1 quasars show AGN-like 
mid-IR colors (Stem et al. 2005; Richards et al. 2009a). The 
positions on the sky of the final samples of QSO-ls and Obs- 
QSOs are shown in Figure EJa), and their distribution in red- 
shift is given in Figure |3] 

The Obs-QSOs are (by definition) optically faint, and so 
few (only 7%) are bright enough to obtain good redshifts from 
MMT optical spectroscopy. AGES targeted objects down to 
a flux limit of / < 20 for sources that are optically extended, 
which is the case for almost all the Obs-QSOs. Therefore 
the vast majority of the Obs-QSO sample has only photomet- 
ric estimates of redshift, derived using an artificial neural net 



technique (iBrodwin et al.ll2006h . Uncertainties on photo-zs 
using this technique for optically-bright quasars are typically 
a, = 0.12(1 +z). However the errors are more difficult to es- 
timate for optically-faint Obs-QSOs, for which there are few 
spectroscopic redshifts for comparison. Photo-z uncertainties 
were discussed at length by H07, with the conclusion that typ- 
ical uncertainties are at most a- = 0.25(1 +z) and are likely 
smaller. Figure|4]shows the photo-zs and spec-zs for the hand- 
ful of Obs-QSOs with spectroscopic redshifts, as well as those 
for the QSO-ls for comparison. The impact of photo-z errors 
in the present clustering analysis are addressed in detail in 
§ 16.31 As discussed in § 16.31 random errors in the photo-zs 
can only tend to decrease the observed clustering amplitude, 
so we expect the present analysis to provide a robust lower 
limit on the clustering of the Obs-QSOs. 

Since the primary aim of this analysis is to compare the 
clustering of quasars with and without obscuration by dust, it 
is imperative that the samples are otherwise matched in key 
properties such as redshift and luminosity. We show in Fig- 
ure |3] that the redshift distributions of the two samples are 
similar, and we obtain bolometric luminosities (Lboi) for the 
quasars by scaling from the rest-frame 8 /im luminosity. We 
compute the flux at rest-frame 8 /im by extrapolating between 
the fluxes at 8 and 24 jiva in the observed frame, and use this 
flux to obtain the monochromatic luminosity vLi, at 8 /im. 
We then mul tiply by a lumi n osity- dependent bolometric cor- 
rection from iHopkins et al.l (120071) . which ranges from fac- 
tors of «8 to 1 1, in order to obtain Lboi- Visual inspection of 
the Spitzer data shows that essentially all of the quasars have 
broadly power-law SEDs at these wavelengths, indicating the 
rest-frame 8 /im emission is indeed dominated by the AGN. 
We note that 49 quasars lie outside the region covered the 
MIPS 24 /im observations, while 26 (~3%) of those inside 
the MIPS area are not detected at 24 fim. For these 75 ob- 
jects, we use the estimates of Lboi derived from the rest-frame 
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Fig. 3. — Redshift distributions for the photometric galaxy sample between 
at 0.5 < z <2 (green dotted line), and the unobscured (blue soUd line) and 
obscured (red dashed line) quasar samples of H07. The histogram for galaxies 
is multiplied by 1/300 so that the distribution can be directly compared to 
that of the AGNs. The redshift range 0.7 < z < 1.8 for which the correlation 
analysis is peiformed is shown by the shaded area. Redshift estimates for the 
galaxies and most of the Obs-Q SOs are derive d from the photometric redshift 
calculations using the method of iBrodwin et aUt20 06). The QSO-ls included 
in the correlation analysis have spectroscopic redshifts from MMT/Hectospec 
observations. 



2 /im luminosity as in § 4.6 of lHOlR . 

The distributions in Lboi are almost identical for the QSO- 
1 and Obs-QSO samples, as shown in the top panel of Fig- 
ure [Tib). The median and dispersion in logLboi (erg s"') is 
(45.86, 0.37) and (45.83, 0.39) for QSO-ls and Obs-QSOs, 
respectively, indicating that the two samples are very well 
matched in bolometric luminosity. For completeness, we note 
that if we use the Lboi estimates derived from rest-frame 2 /im 
in H07 and restrict our analysis to QSO-1 and Obs-QSO sam- 
ples that are matched in Lboh this has a negligible effect on the 
clustering results. 

3.2. Galaxy sample 

The sample of 256,124 galaxies is selected from the deeper 
SDWFS IRAC observations, with a flux Umit [4.5] < 18.6. 
The galaxies are selected to have best estimates of photomet- 
ric redshift between 0.5 and 2, with an average photo-z of 
(z) = 1.09. The sample includes an optical magnitude cut of 
/ < 24 to restrict it to optical fluxes for which the photo-zs are 
well-calibrated. To eliminate powerful AGNs, we have also 
exclu ded any object dete cted in 5 ks Chandra X-ray observa- 
tions (iKenter et al.ll2005h or with 5a SDWFS detections in all 
four IRAC bands and colors in the lStern et"an (12005 ) AGN se- 
lection region. The exclusion of AGNs from the galaxy sam- 
ple removes only 6,979 objects and has negligible effect on 
the results. 

The distribution on the sky of the 256,124 galaxies are 
shown in Figure |2b), and their distribution in photometric 
redshift is shown in Figure |3] Photome tric redshifts a re ob- 
tained using an updated version of the Brod win et al.l (120061) 
algorithm, which is based on template fitting to the optical-IR 

In general, the Li,^] derived from the rest-frame 2 ^m luminosity as 
used in H07 (which did not make use of the 24 /jm data) broadly matches 
that obtain from the extrapolated 8 /.tm flux. However, the median Ltoi ob- 
tained from 2 /.tm is smaller for the Obs-QSOs than for the QSO- 1 s by RiO. 1 5 
dex, primarily because the Obs-QSOs have somewhat redder mid-lR SEDs 
consi stent with the nuclear emission being reddened by dust (e.g.. lHaas et all 
12003) . 
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Fig. 4. — Photometric versus spectroscopic redshift for the QSO-1 (blue 
stars) and Obs-QSO (red circles) IR-selected quasar samples. By selection, 
all the QSO-ls have spectroscopic redshifts, while the Obs-QSOs were gener- 
ally too faint for optical spectroscopy. Only 19 Obs-QSOs have good optical 
redshifts; they are plotted here. 

SEDs. The SED fitting produces a redshift probability den- 
sity function (PDF) for each object, where Piz) represents the 
probability that the object lies at redshift z- (Note that the 
neural net used for the quasar photo-zs does not produce an 
equivalent estimate of the PDF. Thus for the quasars we use 
the best value for the redshift, as discussed in §|4]) P{z) is nor- 
malized such that JP(z)dz = 1. For most galaxies the PDF is 
roughly Gaussian in shape, although often with a broader tail 
toward higher redshift. The typical redshift uncertainties are 
a, ^ 0.1(1 +z), and only a small fraction (ftiO.6%) of galaxies 
show multiple significant peaks in the PDF at different red- 
shifts. Typical galaxy PDFs are shown in Figure|5] 

In addition to the observed galaxy catalog, the correlation 
analysis requires a reference sample of objects with random 
sky positions, in order to compare the observed quasar-galaxy 
pair counts with the number expected for an uncorrected dis- 
tribution. We use a catalog of 8 x 10^ random "galaxies" that 
are assigned to random positions in the regions of good pho- 
tometry, reflecting the spatial selection function for the SD- 
WFS galaxies. 

4. CORRELATION ANALYSIS 

In this section we outline our methods for measuring the 
spatial cross-correlation between quasars and galaxies, the au- 
tocorrelation of the galaxies, and the absolute bias and char- 
acteristic dark matter halo masses. 

4. 1 . Projected correlation function 

To measure the spatial clustering of quasars, we can in prin- 
ciple derive the autocorrelation of the quasars themselves, 
or measure their croii-correlation with a sample of other 
objects (specifically, normal galaxies) at the same redshifts. 
Our quasar sample is too small to obtain sufficiently good 
measurements of their autocorrelation function. However, 
cross-correlation with galaxies (of which there are ~ 300 
times as many objects in the Bootes data set) allows far 
greater statistical power. Further, cross-correlation requires 
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Fig. 5. — Sample probability density functions for three galaxies in the SD- 
WFS sample. Dotted lines show the "best" (peak) redshifts for each galaxy. 
The redshift of a sample quasar is shown by the dashed black line. Note that 
for the two lower-redshift galaxies, the radial distance between the "peak" 
redshift of the galaxy and the quasar redshift are far too large for them to 
be physically associated. However, because of the uncertainty in the galaxy 
redshifts (shown by the PDFs), there is a non-negligible probability that the 
galaxies lie close to the radial distance of the quasar. 



knowledge only of the selection function for the galaxies, 
which is generally better constrained than that for AGNs. 
Cross-correlations of AGNs with galaxies have proved an 
effectiv e techni que in a number of previous s tudies (e.g. , 
Croom et al., ,2004t lAdelberger & SteideH 120051: iSerber et al 
2006t ILietalJ 120061: ICoiletaLI 120071: IWake et alJ 1200: 



Coil et all 12009^, 'Hickox 'etan 



2009t iMan delbaum et al 



2009r 



Padmanabhan 
20091: iMounti-ichas et al. 200' 



et alJ 



Donoso et al.li2010l:lKrumpe et al.|[201Q) . 

For the present analysis, the uncertainties in the galaxy 
photo-zs restrict our ability to perform a full three- 
dimensional clustering analysis. However, making use of the 
quasar redshifts and the galaxy photo-z information, we can 
derive a projected spatial correlation function (wp(R), with R 
in comoving Mpc) that has both higher signal-to-noise, 
and a more straightforward physical interpretation than, for 
example, the purely angular correlation function io{6). 

The two-point correlation function ^(r) is defined as the 
probability above Poisson of finding a galaxy in a volume el- 
ement dV at a physical separation r from another randomly 
chosen galaxy, such that 



dP = n[l+^(r)]dV, 



(1) 



where n is the mean space density of the galaxies in the sam- 
ple. The projected correlation function Wp(R) is defined as the 
integral of ^(r) along the line of sight. 



Jo 



(2) 



where R and tt are the projected comoving separations be- 
tween galaxies in the directions perpendicular and parallel, 
respectively, to the mean line of sight from the observer to the 
two galaxies. By integrating along the line of sight, we elim- 
inate redshift-space distortions owing to the peculiar motions 
of galaxies, which distort the line-of-sight distances measured 
from redshifts. Wp(R) has been used to measu re correla- 
tions in a number of surveys, for example SDSS (Zehavi et al.' 
12005; .Li et al...2006; .Mvers et al...2009; .Krumpe et al. .2010.) . 



2SLAQ (IWakeet al.ll2008h. DEEP2 dCoil etalj|2007l |2QQ1 
120091) . Bo otes ( Hickoxetal. 2009; Starikova et al. 20m 
COSMOS d GiliTet al.«2009 ) and GOODS (.Gilh et al.n2007tl) . 

In the range of separations 0.3 ^ r <5Q Mpc, ^(r) for 
galaxies and quasars is roughly observed to be a power-law. 



ar) = (r/ro)- 



(3) 



For sufficiently large Tr^ax such that we average over all line- 
of-sight peculiar velocities, Wp{R) can be directly related to 
^(r) (for a power law parameterization) by 



Wp{R) = R 



7 r(i/2)r[(7-i)/2] 
r(7/2) 



(4) 



We use Equation ^ to obtain power-law parameters for 
the observed correlation functions, to facilitate straightfor- 
ward comparisons to other works. However, we note that a 
number of recent studies have shown evidence for separate 
terms in the correlation function owing to pairs of galaxies 
found within a single dark matter halo (the "one-halo" term), 
and from pairs in which each galaxy is in a different halo^the 
"two-halo" terrn ; e.g. , Zehavi et al. 2004; Zheng et al. (20071 
ICoil et al.l l2008t iBrown et al.. .2008,; Zheng et al. ,20091). A 
halo occupation distribution (HOD) analysis accounting for 
both the one- and two-halo terms can provide valuable con- 
straints on the distribution of objects within their dark matter 
halos, however a full HOD calculation is beyond the scope of 
the present analysis. 

To measure Wp{R) for the quasar-galaxy cross-correlation, 
we employ the method developed by iM09 . This technique 
makes use of the full photo-z PDF for every galaxy, to weight 
quasar-galaxy pairs based on the probability of their being as- 
sociated in redshift space. We de scribe the formalism briefly 
here, and refer the reader to lM09l for further details. 

4.2. Cross-correlation method 

For a set of spectroscopic quasars all at the same comoving 
distance x* from the observer, the angular cross-correlation 
between the (spectroscopic) quasars and (photometric) galax- 
ies can be expressed in terms of the physical transverse co- 
moving distance by (e.g.. Shanks et al. 1983); 



weiR) ■■ 



Nr DqDc(R) 
Ng DqRc(R) 



1, 



(5) 



where R is the projected comoving distance for a given an- 
gular separation 6, such that R = x*0. Nq and Nr are the 
total numbers of photometric galaxies and random galaxies, 
respectively, and DqDc and DqRq are the number of quasar- 
galaxy and quasar-random pairs in each bin of R. 

Defining the radial distribution function for the full galaxy 
sample as fix), where J f{x)dx = 1- and assuming that f(x) 
varies slowly at the redshifts of interest, then the angular cor- 
relation function wg{R) is related to the projected real space 
correlation function Wp{R) by 



we{R) = f{x*)wp{R) 



(6) 



(for a derivation see § 3.2 of iPadmanabhan et al.l 120091) . As 
discussed in detail in M(j§, we can generalize the analysis 
such that the contribution to WpiR) is calculated individually 
for each quasar-galaxy pair, with fij defined as the average 
value of the radial PDF /(x) for each photometric object /, in 
a window of size Ax around the comoving distance to each 
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spectroscopic object j. We use Ax = 100 /;"' Mpc to effec- 
tively eliminate redshift space distortions, although the results 
are insensitive to the details of this choice. 

In thi s cas e of weighting by pairs, we obtain, as in Equation 
(13) oflMQSt 



w 



/nX AT A7 \ ^ DqDc(R) \ ^ 



'J 



where 



(8) 



We refer the reader to § 2 of lM09l for a detailed derivation and 
discussion of these equations. 

M09 use Equation (|7]) to compute the cross-correlation be- 
tween spectroscopic and photometric quasars from SDSS in 
the relatively narrow redshift bin 1 .8 < z < 2.2, correspond- 
ing to comoving distances 3400 < x < 3800 /i"' Mpc; they 
obtain a cross-correlation le ngth rp = 4.56 ± 0.48 /i"' Mpc, 
assuming 7 = 1.5. (Note that lM09l do not derive 7 = 1.5, they 
assume it purely to describe their method. Higher values of 
gamma are typically obtained in the recent literature, and we 
obtain 7 « 1 .8 in the present work.) Our quasar sample spans 
a comparatively larger range in redshift (0.7 < z < 1.8, corre- 
sponding to 1750 < X < 3400 /z"' Mpc). 

We evaluate Equation by calculating the DqDc /DqRc 
term individually for each quasar. That is, for each quasar and 
each bin in separation R, we sum the redshift weights c, ^ for 
galaxies in the given range of distance from the quasar, and 
divide by the number of random galaxies in the same distance 
range (note that this implies Nq = 1 in Equation|7]i. The advan- 
tage of this procedure is that it consists of a simple sum and 
accounts exactly for the comoving distance to each quasar. 
However, the calculation is limited by shot noise on small 
scales where we have small numbers of quasar-galaxy and 
quasar-random pairs. To check that this does not significantly 
affect the results, we also divide the quasar sample into bins of 
width Az = 0.1 (over which the comoving distance variations 
are small enough that there is little mixing between bins in R), 
and calculate the DqDc/DqRq term for all the Nq quasars in 
the bin. We then average the Wp(R) values for the different 
bins to obtain a mean Wp{R) over the redshift range of inter- 
est. The resulting clustering amplitude differs by <10% (and 
in the majority of cases, < a few percent) compared to evalu- 
ating Equation (Q treating each quasar separately. The choice 
of method does not affect any of our conclusions, but to ac- 
count for these differences we conservatively include an addi- 
tional 10% systematic uncertainty on the measurement of the 
clustering amplitude. Finally, we emphasize that we are av- 
eraging Wp{R) over the whole redshift range of 0.7 < z < 1.8. 
The validity of this procedure depends on the fact that the ob- 
served Wp{R) varies slowly in the redshift range of interest, 
which we verify explicitly in § 16.21 

4.3. Galaxy autocorrelation 

To estimate dark matter halo masses for the quasars, we 
calculate the relative bias between quasars and galaxies from 
which we derive the absolute bias of the quasars relative to 
dark matter. As discussed below, calculation of absolute bias 
(and thus halo mass) requires a measurement of the autocor- 
relation function of the SDWFS galaxies. The large sample 
size enables us to derive the clustering of the galaxies accu- 
rately from the angular autocorrelation function oj{9) alone. 
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Fig. 6.— Redshift distribution for the full sample of 256,124 SDWFS 
galaxies (dotted line) and the sample of 151,256 galax ies s elected to match 
the redshift overlap with the quasars, as described in j} |4.3l (solid line). The 
second gal axy sample is used to derive the angular autocorrelation of the 
galaxies (§ I4.3t , as well as their angular cross-correlation with the quasars 

(§113). 



Although we expect the photometric redshifts for the SDWFS 
galaxies to be well-constrained (as discussed in § I3.2| |. by us- 
ing the angular correlation function we minimize any uncer- 
tainties relating to individual galaxy photo-zs for this part of 
the analysis. The resulting clustering measured for the galax- 
ies has much smaller uncertainties than that for the quasar- 
galaxy cross-correlation. To save computation time, for the 
galaxy autocorrelation analysis we use a significantly smaller 
random catalog with only 5 x 10^ random "galaxies". This 
likely introduces some additional shot noise into the calcula- 
tion of oj(0), however since the resulting uncertainties are still 
far smaller than those for the quasar-galaxy cross-correlation, 
they are more than sufficient for the present analysis. 

We calculate the angular au tocorrelation function uj(d) us- 
ing the lLandy & Szalayl (Il993h estimator: 



uj(e)= —{DD-2DR + RR), 
RR 



(9) 



where DD, DR, and RR are the number of data-data, data- 
random, and random-random galaxy pairs, respectively, at a 
separation 6, where each term is scaled according to the total 
numbers of quasars, galaxies, and randoms. 

The galaxy autocorrelation varies with redshift, owing to 
the evolution of large scale structure, and because the use of a 
flux-limited sample means we select more luminous galaxies 
at higher z- This will affect the measurements of relative bias 
between quasars and galaxies, since the redshift distribution 
of the quasars peaks at higher z than that for the galaxies and 
so relatively higher-z galaxies dominate the cross-correlation 
signal. To account for this in our measurement of galaxy auto- 
correlation, we randomly select galaxies based on the overlap 
of the PDFs wi th th e quasars in comoving distance (in the 
formalism of § 14.21 this is fij for each galaxy, averaged all 
quasars). We select the galaxies so their distribution in red- 
shift is equivalent to the weighted distribution for all galaxies 
(weighted by (fij))- The redshift distribution of this galaxy 
sample is shown in Figure|6] We use this smaller galaxy sam- 
ple to calculate the angular autocorrelation of SDWFS galax- 
ies. 
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4.4. The integral constraint 

In fields of finite size, estimators of the correlation func- 
tion based on pair counts are subject to the integral constraint, 
which can be expressed as (.Groth & Peebles,. 1977b 



(10) 



where 6' 12 is the angle between the solid angle elements 
and dil2 and the integrals are over the survey area. If the 
number density fluctuations in the volume are small, and the 
angular correlations are smaller than the variance within the 
volume, then to first order the correlation function is simply 
biased low by a constant equal to the fractional variance of the 
number counts. A straightforward way to remove this bias is 
to add to the observed ljj{6) the term 



= ^ / / w(6'i2)df7idr22, 



(11) 



where VI is the area of the survey region. The value of 
n^ujn, where n is an estimate of the mean number of galax- 
ies per unit area, is the contribution of clustering to the vari- 
ance of the galaxy nu mber counts dGroth & PeeblesI 1 19771: 
lEfstathiou et al.l 1199 Ih . Evaluating Equation dTTT l for the 
Bootes survey area and the typical slope of the uj{9) for the 
objects considered here, we obtain ojq rs 0.03w(1'). We es- 
timate a;(l') by interpolating the observed 10(6), then add ujn 
to uj(9) before performing model fits. For the projected real- 
space correlation functions Wp(R) (which is ultimately derived 
from individual estimates of coiO), as in Equation |5]l, we per- 
form an approximate correction for the integral constraint . We 
determine the value of Wp{R) at the physical scale (typically 
0.5-1 Mpc) corresponding to 1' for each quasar, and add 
the average of these estimates (multiplied by 0.03) to the ob- 
served Wp{R). These corrections increase the observed clus- 
tering amplitude by «10%, but have little effect on our overall 
conclusions. 

4.5. Uncertainties and model fits 

Ideally, uncertainties in Wp{R) and 10(6) would be deter- 
mined by calculating the correlation function for various ran- 
dom realizations of mock IR-selected quasar and galaxy sam- 
ples, for example by populating dark matter A^-body simula- 
tions. In the absence of such mock catalog, we instead de- 
termine uncertainties in Wp{R) directly from the data through 
bootstrap resampling. 

In a standard bootstrap analysis, the survey volume is di- 
vided into A^sub subvolumes, and these subvolumes are drawn 
randomly (with replacement) for inclusion in the calculation 
of the correlation function. Owing to the relatively small size 
of the field compared to large surveys such as SDSS or 2dF, 
we are only able to divide the field into a small number of sub- 
volumes (we choose A^sub = 8). The width of one sub volume 
corresponds to Ri50 /i"' Mpc at z = 1 .2, so that correlations be- 
tween the subvolumes should be relatively weak. (We verify 
explicitly that using a larger A^sub = 22 has no significant ef- 
fect on the results.) For each bootstrap sample draw a total of 
3A^sub subvolumes (with replacement), which has been shown 
to best app roximate the intrinsic uncertainties in the clustering 
amplitude (iNorberg et"ani2009i) . We then re-calculate Wp(R) 
including only the subvolumes in the bootstrap sample. For 
the calculations of Wp(R) we use 10,000 bootstrap samples, 
for which the uncertainties at each scale converge to better 



than 1%. (To save computing time, we limit the analysis to 
2000 bootstrap samples for the angular correlation analyses, 
for which the uncertainties converge to within ?al.5%.) 

This bootstrap technique works well for the galaxy auto- 
correlation, for which we have a large number of objects and 
the uncertainty is dominated by the clustering of the sample 
rather than counting statistics. However, for the quasar-galaxy 
cross-correlation the bootstrap analysis results in very small 
errors that are significantly smaller than the observed scatter 
between points. This appears to be caused by the fact that, ow- 
ing to the small quasar samples of only a few hundred objects, 
the uncertainties are dominated by shot noise that is not fully 
characterized by randomly selecting entire subvolumes. To 
account for the shot noise, we therefore take the sets of 3A^sub 
bootstrap subvolumes and randomly draw from them (with re- 
placement) a sample of objects (quasars or galaxies) equal in 
size to the parent sample; only pairs including these objects 
are used in resulting cross-correlation calculation. This pro- 
cedure yields a good estimate of the shot noise (the resulting 
~ 1) while also accounting for covariance due to the large- 
scale structure. 

When fitting power-law models to the observed correlation 
functions, we compute parameters by minimizing x^, taking 
into account covariance between different bins in R. From the 
bootstrap analysis, we can estimate the covariance matrix Cij 
by 



1-N 



k=l 



w''p(Rj)-Wp(Rj) 



(12) 



where Wp(Ri) and Wp{Rj) are the projected correlation func- 
tion derived for the A:-th bootstrap samples, A^ is the total num- 
ber of bootstrap samples, and Wp{R) is the correlation function 
for the full sample. This formalism is equally valid for bins 
of angular separation 6 in calculations of ll!{9). The la un- 
certainty in each bin in R is the square root of the diagonal 
component of this matrix (cr, = (C,,)^/^). 
Taking into account covariance, x^ is defined as 

A'bins Nt,„ 



xC-/{wp(Rj)- 



w 



model 



(13) 



where C,^' is the inverse of the covariance matrix Cij. We de- 
termine best-fit parameters by minimizing x^, and derive Icr 
errors in each parameter by the range for which Ax^ = 1 . As 
a check, we also estimate parameter uncertainties by calcu- 
lating best-fit parameters for each of the bootstrap samples 
and calculating the variance between them; this obtains al- 
most identical estimates of the errors. Further, we note that 
if we use only the diagonal terms in the covariance matrix in 
determining x^, the variation in the best-fit parameters is sig- 
nificantly smaller than the statistical uncertainties, indicating 
that the precise details of the covariance matrix are relatively 
unimportant. 

We also note that while in principle the SDWFS field is 
large enough to enable measurements of clustering up to ~50 
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Mpc atz ^ 1, we limit the analysis to scales < 12 /;"' Mpc, 
because of edge effects that skew the correlation function on 
large scales but have minimal effect on smaller scales. An 
investigation of this effect is given in the Appendix. 

4.6. Power law fits to angular correlation functions 

For the projected real-space quasar-galaxy cross- 
correlation analysis, we fit power-law models of Wp{R) 
using Equation (|4]i. We also fit power laws to the an- 
gular correlation functions (both galaxy autocorrelations 
and quasar-galaxy cross-correlations), using the simple 
expression 

ujie)=Ae-^. (14) 

For meaningful comparison to other clustering measurements 
obtained using samples with different distributions in redshift, 
we wish to convert the observed parameters A and 5 to the 
real-space ro and 7 as defined in Equation (|3]l. Inverting Lim- 
ber's equation, the conversion between these par ameters can 
be computed analytically (here we follow § 4.2 of lMvers et al.l 
I2OO6; for the full derivation see ,Peebles..l98(A : 

7 = 5+1 (15) 

.jj /(r(djVi/dz)(cW2/dz)£,x^-"dz 
''"''^ [S^{m/Az)Az] [/„°°(dA^2/dz)dz] ° ^ 

where = r(0.5)r(0.5[7- l])/r(0.57), T is the gamma 
function, x is the radial comoving distance, dA^i,2/dz are the 
redshift distributions of the samples (for an autocorrelation 
dNi/dz = dA^2/dz), and E, = Hjc = dz/dx- The Hubble pa- 
rameter //, can be found via 

H^ = H^[Q,„(l+z)' + nA]. (17) 

Equation ( fT6b assumes no evolution with redshift in the 
clustering of the sample (equivalent to the implicit assump- 
tion made in fitting Wp(R) with Equation |4]l. For each angu- 
lar correlation analysis, we derive A and S from the observed 
uj(9) and then obtain the corresponding 7 and ro from Equa- 
tions ( fTsT i and ( fT6] l. 

4.7. Absolute bias and dark matter halo mass 

The masses of the dark matter halos in which galaxies 
and quasars reside are reflected in their absolute clustering 
bias relative to the dark matt er distribution. To determine 
absolute bias (following e.g., iMvers et"an 120071: ICoil et afl 
12008; Hickox et al. 2009) we first calculate the two-point au- 
tocorrelation of dark matte r as a funct i on of redshift. We 
use the HALOFIT code of ISmith et"aD (l2003h to determine 
the nonlinear-dimensionless power spectrum A^^C^:^) of the 
dark matter assuming our standard cosmology, and the slope 
of the initial fluctuation power spectrum, F = Vl„,h = 0.21. The 
Fourier transform of the A^l(A;,z) gives us the real-space cor- 
relation function ^(r), which we then integrate to tt = 100 /;"' 
Mpc following Equation (|2]i to obtain the dark matter pro- 
jected correlation function vv™(/?,z). The uncertainty in the 
DM power spectrum obtained from HALOFIT is ^ 5%; this 
corresponds to a systematic uncertainty ^0.05 dex in Mhaio, 
but does not impact the relative halo masses of the different 
subsamples. 

To derive quasar absolute bias from the projected real-space 
correlation function, we average the w^^(R) over the redshift 



distribution of the sample, weighted by the overlap with the 
galaxy PDFs. The overlap of each quasar with the galaxy 
PDFs is given by 

W, = ^/,V (18) 

i 

and the corresponding Wp{R) for the dark matter is given by 
w™(/?) = ^ W,w™(«,z,)/ E ^' (19) 

where z, is the quasar redshift. 

The redshift distributions for the QSO-ls and Obs-QSOs 
are essentially identical (the resulting w™(/?) values for the 
two samples differ by <2% on all scales) so for simplicity 
we use the same w^^(R) (defined for the QSO-ls) for both 
sets of quasars. We obtain the bias by calculating the aver- 
age ratio between the best-fit power-law model and w™(/?) 

over the range of scales of 1-10 Mpc, for which w™(7?) 
corresponds closely to a power law and is dominated by the 
two-halo term. The observed clustering amplitude relative to 
the dark matter corresponds to bgbc, where bg and be are the 
absolute linear biases of the quasars and SDWFS galaxies, re- 
spectively. 

To measure be from the galaxy autocorrelation function (or 
bgbc from the quasar-galaxy angular cross-correlation, de- 
scribed in § 16.1b . we require an estimate of the correspond- 
ing uj{9) of the dark matter. To obtain uj{9) we use Lim- 
ber's equation to proje ct the power spectrurn Al,, ( k,z) into 
the angular correlation (Limber 1953: Peeblesl ll980l iPeacockl 
J991; Baugh & Efstathiou 1993). Specifically, we perform 
a Mon te Carlo integration of Equation (A6) of iMvers et al.l 
( 120071) to obtain uj(9) for the dark matter The key param- 
eter in this equation is (dA^c/dz)^ where dNc/dz is the red- 
shift distribution of the galaxies. We calculate dNc/dz from 
the sum of the PDFs of the galaxies for which we perform 
the autocorrelation. In deriving the dark matter lu{9) for 
the quasar-galaxy cross-correlation, we replace (dN/dzf' with 
(d^Vg / dz)(dA^c / dz) where d^Vg / dz is the distribution of quasar 
redshifts. For each angular correlation analysis we compute 
the average ratio between the best-fit power law model and 
the dark matter bj(9) on scales I'-IO', where ijj(9) is dom- 
inated by the two-halo term. This ratio yields b^ for galaxy 
autocorrelations or bgbc for quasar-galaxy cross-correlations. 

Finally, we use bg and be to estimate the characteristic 
mass of the da rk matter ha los hosting each subset of galax- 
ies or quasars. ISheth et aH (12001) derive a relation between 
dark-matter halo mass and large-scale bias that agrees well 
wi th the results ofcosmological simulations. We use Eqn. (8) 
of ISheth et al.l (120011) to convert bahs to Mhaio for the mean 
redshift of each subset of obj ects. If we use a d ifferent re- 
lation between Zjabs and Mhaio linker et al.ll20 05^. we obtain 
estimates for Mhaio that are similar, although slightly larger by 
0.2-0.3 dex; these differences do not significantly affect our 
conclusions. 

We note that to estimate Mhaio, we have performed fits to 
the observed Wp(R) on scales of 0.3-12 Mpc. In princi- 
ple the dark matter and galaxy correlation functions can have 
somewhat different shapes such that the bias depends on the 
range of scales considered. If we limit the fits on scales 1-12 
/i"' Mpc, the results change by ^5%, but with slightly larger 
uncertainties. We also note that our estimates of Mhaio are rel- 
atively insensitive to our choice of trg . If we change from 
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Fig. 7. — The projected quasar-galaxy cross-con'eladon function (derived using Equation |7) for the QSO-1 (blue stars) and Obs-QSO (red circles) samples. 
Uncertainties are estimated from bootstrap resampling. Data points for the two quasar types are slightly offset for clarity. Power-law fits to Wp(R) are shown as 
solid lines in blue (QSO-ls) and red (Obs-QSOs), and the projected correlation function for dark matter is shown by the dotted gray line. Fits are performed over 
the range in separation of 0.3 < R < 12 A"' Mpc. 



0.84 to 0.8 (as favored by th e more recent WMAP cosmol- 
ogy, e.g. ISpergel et al.ll200'7h our Mhaio estimates for quasars 
and galaxies increase by 0. 1 dex. 

5. RESULTS 

In this section we discuss the results of the correlation anal- 
ysis and the characteristic dark matter halo masses for galax- 
ies and quasars. We first calculate the cross-correlation of 
the full QSO-1 and Obs-QSO samples with SDWFS galaxies. 
The resulting Wp(R) values and best-fit models are shown in 
Figure [T] and fit parameters are given in Table [T] For both 
sets of the quasars the observed real-space projected cross- 
correlation is highly significant on all scales from 0.1-12 
/i"' Mpc, and the power law fits return 7 ft; 1.8, similar to 
many previous correlat i on function meas urements for quasars 
(e.g.. ICoil efal] 120071: iRoss et alj 120091) and galaxies (e.g., 
IZehavi et al.ll2005l: ICoil et alj|2008l) . "fhe best-fit parameters 
are ro = 5.4 ± 0.7 /z"' Mpc, 7 = 1 .8 ± 0. 1 for the QSO-ls, and 
ro = 6.4±0.8/!-' Mpc, 7= 1.7 ±0.1 for the Obs-QSOs. The 
results indicate that the cross-correlation of the Obs-QSOs 
with galaxies is somewhat stronger than that for the QSO-ls. 
The corresponding values of bghc are given in Table [T] 

As a check, we also perform power law fits to Wp(R) but 
leaving the slope fixed to 7 = 1.8, which corresponds to the 
slope of the Wp{R) for the dark matter This also yields a sig- 
nificant difference in the clustering amplitude, although some- 
what smaller, with ro = 5.3 ± 0.6 and ro = 6.0 ± 0.6 /;"' Mpc 
for the QSO-ls and Obs-QSOs, respectively. (Note that the 
formal uncertainties in ro here are smaller than for the above 
results because they do not account for covariance with 7.) 

To obtain the absolute bias of SDWFS galaxies (be) in or- 
der to extract the quasar bias from the cross-correlation 
results, we next derive the autoco rrela tion of SDWFS galax- 
ies for the sample described in § 14.31 The observed u>(9) is 



shown in Figure [8] along with the correlation function for 
dark matter, calculated as discussed in § 14.71 Fit parame- 
ters are given in Table [1] The power-law model fits well on 
the chosen scales of l'-12', although there is a clear excess 
corresponding to the one-halo term at 6* < 1', a s is common 
in gal axy autocorrela tion measurements (e.g., i Quadri et alj 
2008; Kim et alj|201 Ih . The best-fit power law parameters are 
ro = 4.7 ± 0.2 and 7 = 1 .67 ± 0.05, and the ratio of the best-fit 
power law to the dark matter uj(9) yields b^ = 2.79 ± 0.16 or 
be =1.61 ±0.05. 

This accurate value for be allows us to estimate bQ for 
both types of quasars, based on the cross-correlation measure- 
ments. We obtain bg = 2.17 ±0.55 and 3.06 ±0.70, for QSO- 
ls and Obs-QSOs, respectively. Convert ing these to dark mat- 
ter halo masses using the prescription of Sheth et all (1200 ID as 
described in §|4J] we aiTive at log(Mi^^io[hM^]) = 12.7+|];^ 
and 13.3^1 for QSO-ls and Obs-QSOs, respectively. The 
difference is marginally significant (« Icr, although as we dis- 
cuss below, the Obs-QSO clustering may represent only a ro- 
bust lower limit). 

For direct comparison with other studies that directly mea- 
sure the quasar autocorrelation, it is useful to present the 
quasar clustering in terms of effective power law parameters 
for their autocorrelation. Assuming linear bias, the quasar 
autocorrelation can be inferred from the cross-correlation by 
^QQ - ^Qc/^cc (e.g.. lCoil et al.ll2009h . We can therefore use 
the power law fits to the quasar-galaxy cross-correlation and 
galaxy autocorrelation to derive an effective ro and 7 for 
the quasar autocorrelation. This yields ro = 6.1^[ g /;"' Mpc 
and 7 = 1.9!;];^ for the QSO-ls and ro = 8.8!|!^ /z"' Mpc and 
7 = 1.7 ±0.2 for the Obs-QSOs. The autocorrelation ampli- 
tude and Mhaio for QSO-ls are in excellent agreement with 
previous estimates for unobscured quasars, while the best-fit 
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amplitude for Obs-QSOs is higher than most previous mea- 
surements of quasar clustering. We compare these results to 
previous work and discuss possible interpretations in § |7] 

6. VERIFICATION 

In this section we perform several tests to verify the validity 
of the clustering analysis outlined in §2] We first calculate the 
quasar-galaxy cross-correlation using a simple angular corre- 
lation function, minimizing dependence on the photometric 
redshifts. We then check for variation in the observed clus- 
tering over the redshift range of interest and confirm that any 
variation is relatively weak. Finally, we estimate the effects 
of uncertainties on the photometric redshifts on the observed 
real-space clustering amplitude for the Obs-QSOs. These 
checks confirm that our projected correlation analysis pro- 
vides a robust estimate of the quasar-galaxy cross-correlation. 

6.1. Angular cross-correlation 

We first calculate the cross-correlation of quasars and SD- 
WFS galaxies using a simple angular clustering analysis, and 
check whether the corresponding absolute bias is consistent 
with that derived from the more sophisticated Wp(R) calcula- 
tion. To calculate the 10(6) we use an estimator corresponding 
to Equation (O but for cross-correlations: 



1 

RR 



(DqDg-DqR-DgR+RR), 



(20) 



where each term is scaled according to the total numbers of 
galaxies and randoms. To maximize the signal-to-noise ratio 
by cross-correlating objects associated in redshift space, the 
galaxies include only the redshift-matched SDWFS sample of 
151,256 objects described in 14. 31 Uncertainti es ar e estimated 
using bootstrap resampling as described in § 14.51 We fit the 
observed cross-correlation with a a power law as described 
in § 14.61 Owing to the limited statistics which provide only 
very weak constraints on the power law slope, we fix (5 = 0.8 
(corresponding to real-space 7 = 1.8). 

The resulting cross-correlations and scaled dark matter fits 
are shown in Figure |9] and fit parameters are given in Ta- 
ble [U The estimates of ro and bgbc are in broad agreement 
between the two estimators, although as may be expected, the 
statistical uncertainties for the angular correlation analysis are 
larger (by ^^50%) than for the real-space analysis with fixed 7. 
Given that the absolute bias derived from the projected corre- 
lation function corresponds broadly to the bias from the nois- 
ier, but simpler angular cross-correlation, we conclude that 
there are no significant systematic effects that skew our esti- 
mate of Wp{R). 

6.2. Variation in Wp{R) with redshift 

Our calculation of the real-space quasar-galaxy correlation 
function over the redshift range 0.7 < z < 1.8 requires that 
Wi,{ R) va ries only slowly between these redshifts, as discussed 
in § 14.21 If the objects reside in similar halos at all redshifts, 
then we may expect ro to change slowly; simulations suggest 
that the typical ro for the autocorrelation of DM halos of mass 
-lO'^-lO'^ /z"' Mq should cha nge by < 0.2 h~^ Mpc between 
z = 0.5 and 2 (see Figure 10 of IStarikova et al.l l2010). To test 
explicitly the redshift variation for the clustering in our sam- 
ple, we re-derived Wp(R) using the method outlined in §|4]but 
selecting quasars over smaller redshift bins of 0.7 < z < 1 .25 
and 1.25 < z < 1.8. Uncertainties are calculated using the 
bootstrap method as for the full quasar samples, and dark mat- 
ter and power-law fits are again performed over the range of 
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Fig. 8. — The angular autocorrelation function of SDWFS galaxies, se- 
lected to match the overlap of the quasars and galaxies in redshift space. Un- 
certainties are estimated from bootstrap resampling. The angular correlation 
function for dark matter, evaluated for the redshift distributions of the galax- 
ies, is shown by the dotted gray line. The power law fit was performed on 
scales l'-12' and is shown as the solid line. The excess in uj{9) aX 9 < I' is 
due to the one-hal o term arising from pairs of galaxie s within the same dark 
matter halos (e.g. JOuadri et al.l2008trKim et alJ20Tll) . 

separations 0.3 < R < 12 h~^ Mpc. We evaluate Wp^(R) as 
in § 14.71 but only including the quasars in the redshift ranges 
of interest. Owing to larger statistical errors and for simple 
comparison to the results over the full redshift range, in the 
power law fits we fix 7 to 1.8. 

The resulting Wp{R) for the separate redshift bins and the 
power law fits are shown in Figure [TO] For the QSO-ls 
we obtain ro = 5.1 ±0.8 /;"' Mpc and 6.3 ±0.7 /i"' Mpc 
for the low- and high-redshift bins, respectively, and for 
the Obs-QSOs we correspondingly obtain ro = 6.2 ±0.7 /i"' 
Mpc and 5.8 ± 1.0 /;"' Mpc. The measured quasar-galaxy 
cross-correlation should be largely independent between the 
two redshift bins. Although the quasars are cross-correlated 
against the same galaxy sample in each bin, the galaxy sam- 
ples will be weighted toward higher and lower redshifts in 
the high- and low-z bins, respectively. For the high and low 
redshift bins, the best-fit ro values bracket those for the full 
redshift samples, and are broadly consistent within the un- 
certainties. Interestingly, the best-fit clustering amplitude for 
the Obs-QSOs increases with redshift while it decreases for 
the Obs-QSOs; however, given the uncertainties we decline 
to speculate on any possible difference in redshift evolution 
between the two subsets. Overall, the results in the different 
redshift bins confirm that any variation in the observed Wp{R) 
is sufficiently weak over the redshift range of interest, so that 
the method outlined in § 14. 2 1 should provide a reasonable esti- 
mate of the average clustering amplitude over the full redshift 
range. 

6.3. Effects of quasar photo-z errors 

The primary uncertainty in our estimate of Wp(R) for the 
Obs-QSOs is the lack of accurate (that is, spectroscopic) red- 
shifts and difficulty in estimating the photo-z uncertainties 
from the neural net calculations. As described in § H] in 
calculating Wp{R) for the Obs-QSO-galaxy cross-correlation, 
we simply assume that Obs-QSOs lie exactly at the best red- 
shifts output by the neural net estimator Any uncertainties 
in the photo-zs or systematic offsets from the true redshifts 
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TABLE 1 
Correlation results 









Power law fit 




Bias"^ 




Halo mass'^ 


Subset 






ro {h~^ Mpc) 7 








(logr'Mo) 



Projected real-space cross-correlation (Wp{R)f 



QSO-1 


445 


1.27 


5.4±0.7 


1.8±0.1 


1.1 


3.63 ±0.92 


2.17±0.55 


12.7+"'' 


Obs-QSO 


361 


1.24 


6 4+0.7 


1.7±0.1 


1.2 


5.11±1.16 


3.06 ±0.70 


13 3+^-^ 








Wp{R) with fixed 7'' 








QSO-1 


445 


1.27 


5.3 ±0.6 


1.8 


1.1 


3.44 ±0.75 


2.06 ±0.45 


12 6+" '' 


QSO-1 (photo-z) 


445 


1.27 


4.9 ±0.6 


1.8 


1.2 


3.05 ±0.67 


1.82 ±0.40 


12.5+"'' 


Obs-QSO 


361 


1.24 


6.0 ±0.7 


1.8 


1.2 


4.38 ±0.84 


2.62±0.51 


1 3 0+^-^ 










Angular correlation {uj(d)f 








galaxies 


151256 


1.10 


4.7 ±0.2 


1.67 ±0.05 


1.1 


2.79±0.16 


1.67 ±0.05 


12.30±0. 


QSO-1 


445 


1.27 


5.6 ±0.8 


1.8 


1.1 


4. 19 ± 1.07 


2.50 ±0.65 


1 3 0+" '' 


Obs-QSO 


361 


1.24 


6.0 ±1.0 


1.8 


1.2 


4.80 ±1.29 


2.87 ±0.77 


13 2+"" 



" Number of objects include in the correlation analysis. For quasar-galaxy cross-correlation, we use the full sample of 256, 124 galaxies (for Wp{R) calculations) 
or 151,256 galaxies (for ui(9) calculations). 

'° Median redshift for the objects included in the correlation analysis. 

Uncertainties in the DM power spectrum introduce an additional systematic eiTor of ~5% in h.jh^ (and corresponding ~0.05 dex in M\^.^\„). Further systematic 
errors in Mhaio of ~0.2 dex are caused by uncertainty in eg and in the conversion from foabs to ^halo, as discussed in § 14.71 However, these do not significantly 
effect the relative halo masses, so these uncertainties is not included here. Note that for fits with fixed 7, uncertainties on ro, bias, and Mhaio do not account for 
covariance with 7 and thus somewhat underestimate the error on the clustering amplitude. 

Real-space projected cross-correlation between quasars and galaxies, calculated as d escribed in §|4] For all Wp{R) calculations, error estimates for ro, bias, 
andMdaio include a 10% systematic uncertainty on the amplitude as described in § 14.21 
° Angular galaxy autocorrelation and quasar-galaxy cross-correlation, calculated as described in §|4] 

jects (the QSO-ls) that have a similar redshift distribution 
and for which the redshifts are known precisely from spec- 
troscopy. We can therefore adjust the redshifts of the QSO-1 s 
and re-calculate Wp{R) to determine how uncertainties or sys- 
tematic shifts affect the observed correlation amplitude. As a 
simple first test, we calculate l\\tWj,{R) for the QSO-ls using 
the photo-zs (as shown in Figure HJ rather than spectroscopic 
redshifts. Figure[TT]shows that the resulting Wp{R) differs lit- 
tle from that obtained using spectroscopic redshifts; the clus- 
tering amplitude for a power-law fit with fixed 7 = 1 .8 is lower 
by 12%. Note that if we allow 7 to float, the average bias for 
the photo-z sample is actually larger by wlO% (owing to a 
slightly flatter slope) but well within the statistical uncertain- 
ties. We conclude that for the QSO-ls, photo-z errors do not 
have a significant effect on our measurements of the cluster- 
ing. 



QSO-1 ★ 




e (arcmin) 



Fig. 9. — The angular' quasar-galaxy cross-coiTelation function for the 
QSO-1 (blue stars) and Obs-QSO (red circles) samples. Uncertainties are 
estimated from bootstrap resampling. Data points for the two quasar types 
are slightly offset for clarity. The angular correlation function for dark mat- 
ter, evaluated for the redshift distributions of the galaxies and the QSO-ls. 
is shown by the dotted gray line, and power law fits (with fixed 5 = 0.8) are 
shown as solid blue and red fines. All fits ai'e performed over the range in 
separation of 1' to 12'. 



could therefore affect the resulting clustering measurement. 
The fact that we obtain very similar estimates of bq for the 
Obs-QSOs from a simple angular cross-correlation analysis 
as from our Wp{R) calculation suggests that uncertainties on 
individual photo-zs do not strongly affect our estimates of the 
quasar bias, as long as the overall distribution in redshifts for 
the Obs-QSOs is correct. However, it is possible that very 
large discrepancies from the true photo-zs, or any systematic 
shift in the redshift distribution, could affect both the esti- 
mate of Wp{R) and the real-space clustering parameters de- 
rived from Luifi). 

To precisely explore the effect of these errors, we take ad- 
vantage of the fact that we have an equivalent sample of ob- 



6.3.1. Random errors 

To explore photo-z errors in more detail, we systematically 
test the effects of Gaussian random errors in the quasar red- 
shifts. For each quasar, we shift the best estimates of redshift 
(spec-zs for QSO-ls and photo-zs for Obs-QSOs) by offsets 
Az/(1 +z) selected from a Gaussian random distribution with 
dispersion cr^/il +z). Using these new redshifts we recalcu- 
late Wp{R), using the full formalism described in § |4] We 
perform the calculation ten times for each of several values 
of (T,/(l +z) from 0.02 up to 0.2 (which smears the redshifts 
across most of the redshift range of interest). To ensure that 
this step does not artificially smear out the redshift distribu- 
tion beyond the range probed by the galaxies, we require that 
the random redshifts lie between 0.7 < z < 1.8; any random 
redshift that lies outside this range is discarded and a new 
redshift is selected from the random distribution. For each 
trial we obtain the relative bias by calculating the mean ratio 
of Wp{R), on scales 1-10 Mpc, relative to the Wp{R) for 
the best estimates of redshift. We then average the ten trials 
at each a, to obtain a relation between relative bias and a,, 
shown in Figure [T2l a). 
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Fig. 10. — Projected quasar-galaxy cross-correlation function in two redshift bins: 0.7 < z < 1.25 and 1.25 < z < 1.8, with symbols shown as in the upper 
light. Data points for the two redshift bins are slightly offset for clarity. Power-law fits with 7 fixed to 1.8 are shown as solid hues, with c olors corresponding to 
the symbols for each redshift bin. Any variation in Wp(R) with redshift is relatively weak, confirming that the method outHned in j) |4.2| should be valid for this 
analysis. Note that for the Obs-QSOs at z > 1.25, we exclude the bin with Wp(R) < 0, which disproportionately affects the fit. (Including this bin decreases the 
clustering amplitude by Ri20%.) Fit parameters are dark matter halo masses for each redshift bin are given in Table[T] 



10000 



QSO-1 (photo-;) * 
QSO-1 ★ 




1.0 

R (h ' Mpc) 
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Fig. 1 1. — Projected quasar-galaxy cross-correlation function for QSO-ls, 
using photometric redshifts (purple open stars) and spectroscopic redshifts 
(blue solid stars). Power-law fits are performed using a fixed 7 = 1.8, and fit 
parameters are dark matter halo masses given in Table[T] Using the less accu- 
rate (photometric) redshifts has relatively small effect on WpiR), decreasing 
the best-fit power law amplitude by f»10%. 

As may be expected, Figure [T2fa) shows that shifting the 
QSO-1 redshifts from their true values causes a decrease in 
the cross-correlation amplitude, as the quasars are preferen- 
tially correlated with galaxies that are not actually associ- 
ated in redshift space. We find a monotonic decrease in rela- 
tive bias with cr-., from « 0.95 for (7r/(l+z) = 0.05 to « 0.8 
for (Tz/Cl +z) = 0.2. Repeating this calculation for the Obs- 
QSOs reveals a very similar trend. The decrease in bias with 
(T,/(l + z) shown in Figure \T2\ a) indicates that such errors 
would affect the measurements of the clustering amplitude by 
at most --20%. 

6.3.2. Systematic shifts 

While the above analysis suggests that random errors in the 
Obs-QSO photo-zs do not strongly affect the observed cluster- 
ing amplitude, it is also possible that systematic uncertainties 
in the photo-z (consistent over- or under-estimates of the red- 



shift) could significantly alter the observed bias. To te st this, 
we shift the redshifts of the quasars as discussed in § 16.3.11 
but in place of random shifts, we compress all redshifts to- 
ward one end or the other of the 0.7 < z < 1.8 range. (This 
procedure allows us to test the effects of systematic shifts in 
redshift while keeping the same overall redshift range.) The 
shift in redshift is defined by a redshift scaling parameter S^, 
such that 



2new ~ 



z+S,(z-Q.7) 
z+S,(1.8-z) 



S,<0 
S,>0 



(21) 



As an additional check we also perform a simple linear off- 
set Az of the redshifts, allowing the redshifts to move outside 
the selection range of 0.7 < z < 1.8. As in § 16.3.11 we use 
these new redshifts to recalculate Wp(R) via the full formalism 
described in §|4] and determine the relative bias on scales 1- 
10 Mpc. Relative bias versus and Az are shown in Fig- 
ures [T2tb) and (c). For the QSO-ls, the peak of the observed 
clustering amplitude is very close to = 0, while shifting the 
redshifts down or up systematically decreases the bias. The 
Obs-QSOs show a similar peak near S, = 0, indicating that 
the Obs-QSO photo-zs are not systematically offset higher or 
lower than the true redshifts by a large factor We note that 
for the Obs-QSOs a slight shift to higher redshifts (Az = 0.1) 
increases the clustering by a small amount («8%). 

Finally, we emphasize that any possible low-redshift con- 
taminants (such as star-forming galaxies, as discussed in § 7 
of H07), will serve only to decrease the observed cluster- 
ing signal, as they will be completely uncorrected in angular 
space with the higher-redshift SDWFS galaxies that lie in en- 
tirely separate large-scale structures. Therefore the observed 
Wp{R) represents a robust lower limit to the clustering ampli- 
tude for the Obs-QSOs. 

7. DISCUSSION 

We have used the IR-selected quasar sample of lH07l to mea- 
sure the clustering amplitude and to estimate characteristic 
dark matter halo masses for roughly equivalent samples of un- 
obscured and obscured quasars. We obtain highly significant 
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Fig. 12. — Tests of the effects of redshift ei Tors on the observed clustering, 
for QSO-ls and Obs-QSOs, as described in § 16.31 (a) Relative clustering am- 
plitude, after shifting the quasar redshifts by an offset drawn from a Gaussian 
random distribution with width ct^. Both types of quasar show a very similar 
monotonic decrease in clusteiing amplitude with c!-/{\ +z), indicating that 
the photo-z estimates for the Obs-QSOs are reasonably accurate, (b) Corre- 
sponding change in bias after scaling the quasar redshifts toward the limits 
of the redshift interval 0.7 < z < 1.8; the scale parameter S- is defined in 
Equation {2\\ . (c) Relative bias after a simple upward or downward shift of 
the quasar redshifts. In both panels (b) and (c), both types of quasar show 
clustering amplitudes that peak near 5; = or Az = indicating that there is 
no large systematic error in the photo-z estimates for the Obs-QSOs. 



detections of the clustering for both samples, with marginally 
stronger clustering for the Obs-QSOs. In this section, we 
compare our results for QSO-ls to previous results on unob- 
scured quasar clustering, speculate on physical explanations 
for possible stronger clustering for Obs-QSOs, and discuss 
prospects for future studies with the next generation of obser- 
vatories. 

7.1. Comparison with previous results 

We compare our observed absolute bias for the mid-IR se- 
lected quasars to the clustering of optically-selected (Type 
1) quasars, which has been well established by a number 
of works. Among the most precise measurements to date 
are studies that have used data from the 2dF and SDSS sur- 
veys. Recently, IRoss et al.l (|2009) measured the evolution 
of the quasar 3-D autocorrelation function based on spectro- 
scopic quasars from SDSS Data Release 5, and compared to 
previous results from spectroscopic samples from the 2QZ 
(ICroom et al.1 120051) and 2SLAQ (da Angela et al. 2008), as 
well as the clustering of pho tometrically-selected quasars 
from SDSS ( iMvers et al.ll2()06h . Figure [Ha) shows the red- 
shift evolution of Unear bias f or Type 1 quasars from these 
studies, taken from Figure 12 of lRoss et al.l (1 2009*). Where ap- 
propriate, we have converted the bias to our adopted cosmol- 
ogy us ing the formalism in the Appendix of Starikova et al] 
(120101) . assuming a shape factor s = \ for simplicity. Fig- 
ure [T3lb) shows the corresponding estimates of characteris- 
tic Mhain derived fr om the linear bias using the prescription of 
iSheth et aTl ( 120011) . The linear bias and halo masses for the 
QSO-ls and Obs-QSOs are shown for comparison. 

It is readily apparent from Figure [T3}a) that the observed 
bias of Type 1 quasars increases with redshift, as discussed in 
§ [U (These results are also consistent with a number of other 
qu asar clustering studies using other data see e.g.. Figure 15 
of lHopkins et"al] 12008 .) The dashed curves in Figure [T3l a) 
show the increase in bias with redshift for halos of constant 
mass, clearly showing that at all redshifts the QSO-ls reside 
in dark matter halos of roughly a few xlO'^ /;"' Mq. Our 
measurement for the QSO-ls is in excellent agreement with 
the evolution in linear bias and roughly constant Mhaio ob- 
served in previous measurements of Type 1 quasar clustering . 

For the Obs-QSOs, the best-fit bias is marginally larger 
(« Id), corresponding to a factor of roughly four difference 
in the characteristic Mhaio between the QSO-ls and Obs-QSOs 
(FigurefTjb). As discussed above, random errors in the photo- 
zs can only decrease the observed clustering amplitude and 
the inferred Mhaio- Thus our measurement of Obs-QSO clus- 
tering represents a lower limit, and it is possible that the true 
Obs-QSO bias is somewhat higher (although the results of §|6] 
suggest the true bias may be larger by at most ^20%). Based 
on this analysis we can make the robust conclusion that the 
Obs-QSOs are at least as strongly clustered as their QSO-1 
counterparts. 

7.2. Physical implications 

Stronger clustering for obscured quasars would have signif- 
icant implications for physical models of the obscured quasar 
population. In terms of unified models, a difference in clus- 
tering between obscured and unobscured quasars would rule 
out the simplest picture in which obscuration is purely an 
orientation effect, but may be consistent with more compli- 
cated scenarios where the effective covering fraction changes 
with environment. Alternatively, if obscuration is caused by 
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Fig. 13. — Comparison of (a) linear bias and (b) inferred dark matter halo mass to previous measurements of Type 1 quasar clustering. Lines in the left plot 
correspond to the linear bias for dark matter halos with logMhaio[/' 'Mq] = 12.0, 12.5, and 13.0 (from bottom to top). The gray points show results from quasar- 
quasar correlation measurements using spectroscopic samples from SDSS (Ross et al.l'2O09t), 2QZ jCroom et al.ll2005h . and 2SLAQ Ida Angela et alj|2008l) . as 
well as the clustering of photometrically-selected quasars from SDSS (Mvers et al. 2006). The blue star and red circle show our measurements for QSO-ls and 
Obs-QSOs, respectively. Our QSO-1 measurement agrees closely with previous work on Type 1 quasar-quasar correlations, while the Obs-QSOs show marginally 
Si 1(T stronger clusteiing, corresponding to a factor of roughly four in M^aio. 



large (~kpc scale) structures, then the processes that drive 
these asymmetries (e.g., mergers, disk instabilities, accretion 
of cold gas from the surrounding halo) may be more com- 
mon in halos of larger mass. Indeed, given that some fraction 
of quasars might naturally be expected to be obscured by a 
"unified"-model torus, any observed differences in clustering 
may reflect even stronger intrinsic dependence of large-scale 
obscuration and environment. 

An intriguing scenario for obscured quasars is that they 
represent an early evolutionary phase of rapid black hole 
growth before a "blowout" of the obscuring material from 
the central regions of the galaxy and the emergence of an 
unobscured quasar (e.g.. Figure 1 of iHopkins etaH l2008h . 
Quasars tend to radiate at large fraction of the Edding- 
ton rate (McLure & DunloD 2004; KollmeieretaD 120061: 
although see Kelly et al. 2010), so that the similar Lboi 
for QSO-ls and Obs-QSOs would imply that they host 
black holes of similar masses. A ny correlation (e.g., 
Ferrarese"2002^, 'Booth & Sch"ayell2010h . even if indirect (e.g., 
Kormendy & Bender 201 ij, between the final masses of black 
holes and those of their host halos would thus suggest that 
our obscured and unobscured quasars would have the same 
^haio, as long as their black holes are near their final masses. 
However, if obscured quasars are in an earlier phase of rapid 
growth and so are in the pr ocess of "catching up" to their 
final mass (e.g.. iKiii^ 1201 0*). then they would have a larger 
A^haio compared to unobscured quasars with the same Mbh- 
In light of recent debate as to whether bl ack holes gener- 



ally grow befor e or after their hosts (e . g.. Peng et al. 



[Alexander et all 2008; W oo et al.ll200l iDecarh et all 



this scenario would imply that black hole growth lags behind 
that of the host halo. 

In any physical picture, a significant difference in cluster- 
ing between obscured and unobscured quasars would also im- 
ply a difference in accretion duty cycles (or equivalently, life- 
time). QSO-ls and Obs-QSOs are found in roughly equal 
numbers, but th e abundance of dark matter halos drops rapidly 
with mass (e.g.. lJenkins et al.ll2001h . thus implying that if one 
type of quasars are found in larger halos then they must be 



longer-lived. With our current results, we are able to rule 
out any model in which obscured quasars are substantially 
less strongly clustered or have shorter lifetimes than their 
unobscured counterparts. With more accurate future mea- 
surements, detailed studies of halo masses and lifetimes for 
obscured and unobscured quasars could place powerful con- 
straints on evolutionary scenarios such as those described 
above. 

7.3. Future prospects 

Our results demonstrate the potential for studying the clus- 
tering of obscured quasars in extragalactic multiwavelength 
surveys, and the marginally significant difference in cluster- 
ing we observe for obscured and unobscured quasars pro- 
vides strong motivation for more precise measurements in 
the future. The two main avenues for progress are improve- 
ments in redshift accuracy and selection of larger samples 
for better statistical accuracy. Upcoming sensitive, wide- 
field multi-object spectrographs will enable efficient measure- 
ments of redshift for large numbers of optically-faint sources 
and so improve calibrations of obscured quasar photo-zs, 
or with large enough samples, enable fully 3-D clustering 
studies. In addition, we will soon have the capability to 
detect many thousands of obscured quasars based on very 
wide-field observations in the mid-IR with the Wide-Field 
Infrared Survey Explorer (Wri ght! r2008h and in X-rays with 
eROSITA (Predehl et al. 2007) or the Wide-Field X-ray Tele- 
scope (iMurray et al.ll2010i) . These data sets will allow us to 
measure obscured quasar clustering with statistical precision 
that is comparable to current measurements of unobscured 
quasars. 

8. SUMMARY 

We have used data from the Bootes wide-field multi- 
wavelength survey to measure the two-point spatial cross- 
correlation between unobscured (QSO-1) and obscured (Obs- 
QSO) mid-IR selected quasars in the redshift range 0.7 < z < 
1.8. The QSO-ls exhibit clustering corresponding to a typ- 
ical Mhaio ~ 5 X lO'^ Mq, similar to previous studies of 
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Fig. 14. — Observed projected cross-correlation between quasars and galaxies derived using Eqn. (blue stars and red circles for QSO-ls and Obs-QSOs, 
respectively), compared to the same quantity but derived for samples of quasars with the redshifts equal to those of the QSO-ls but with randomized sky positions. 
The gray shaded region shows the dispersion between ten different random quasar samples. Note that the congelation between random quasai's and galaxies is 
small relative to the real-quasar galaxy correlation, except on scales > 10-20 ft"' Mpc where both quantities show a hump and then become negative. This feature 
is likely due to edge effects in the finite SDWFS field, as discussed in the Appendix. Because of this artifact we restrict our congelation analyses to R < 12 /j"' 
Mpc. 



optically-selected quasar clustering. We robustly determine 
that the Obs-QSOs are clustered at least as strongly as the 
QSO-ls, with a marginally stronger signal corresponding to 
host halos of mass '--^ 2 x 10'^ M©; the true clustering am- 
plitude could be up to ^-^20% larger owing to photo-z uncer- 
tainties for the Obs-QSOs that can decrease the observed cor- 
relation amplitude. Our results motivate more accurate mea- 
surements of obscured quasar clustering with larger quasar 
samples and more accurate redshifts. If future studies con- 
firm that obscured quasars are more strongly clustered than 
their their unobscured counterparts, this would rule out the 
simplest "unified" models and may provide evidence for sce- 
narios in which rapid obscured accretion represents an evolu- 
tionary phase in the growth of galaxies and their central black 
holes. 
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APPENDIX 

EDGE EFFECTS AND CORRELATIONS ON LARGE SCALES 

The large area of the Bootes survey allows us to measure galaxies at relatively large physical separations; 1 degree corresponds 
to 50 Mpc at z = 1.5. However, when we use Equation (|7]l to calculate the projected real-space correlation function on 
large scales, we find that the Wp{R) flattens on scales R ^ 10-20 /;"' Mpc, corresponding to tens of arcmin, and then becomes 
negative at ^ 1 00 Mpc (Figure fT4li. This behavior is not observed for quasar clustering from other, wider-field surveys 
(e.g. Groom et al]l2005HMyers et al., 2009) . for which the correlation function continues to decrease on larger scales. While the 
integral constraint require that the correlation fun ction becomes negativ e on some scales, in galaxy auto-correlation surveys this 
generally only happens at 7? > 200 Mpc ('e.g..'Eisenstei n et"ani2005h . 

One possibility is that the observed behavior is due edge effects arising from the finite geome try of the SDWFS field , which 
are not taken into account by the simple ^(/?) = DD/DR- 1 estimator in the IM09I formalism (e.g.. Land v & Szalavlll993l) . To test 
this possibility, we re-performed the correlation analysis described in §|4] after randomizing the positions of the quasars on the 
sky within the area of good SDWFS photometry. We performed 10 separate random trials, for which the cross-correlations are 
shown in Figure [14] along with the Wp{R) values for the QSO-ls and Obs-QSOs. It is clear from Figure [14] that on scales <10 
Mpc, the projected cross-correlation between the random quasars and galaxies is small compared to the Wp{R) for the real 
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quasar sample. However, on scales >20 h ' Mpc, both the real and random samples show an increase in Wp{R) which eventually 
becomes negative around R=1QQ Mpc. The quantities of interest in this paper (i.e. absolute bias and dark matter halo mass) 
can be measured by studying the correlations on scales < 12 /i"' Mpc, where the artifacts are small and have negligible impact on 
the fits to the correlation function. For this paper we therefore limit the correlation analyses to those scales. 
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